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We report on devices based on graphene charge detectors (CDs) capacitively coupled to graphene 
and carbon nanotnbe quantum dots (QDs). We focus on back action effects of the CD on the 
probed QD. A strong influence of the bias voltage applied to the CD on the current through the 
QD is observed. Depending on the charge state of the QD the current through the QD can either 
strongly increase or completely reverse as a response to the applied voltage on the CD. To describe 
the observed behavior we employ two simple models based on single electron transport in QDs 
with asymmetrically broadened energy distributions of the source and the drain leads. The models 
successfully explain the back action effects. The extracted distribution broadening shows a linear 
dependency on the bias voltage applied to the CD. We discuss possible mechanisms mediating the 
energy transfer between the CD and QD and give an explanation for the origin of the observed 
asymmetry. 
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I. INTRODUCTION 

Charge sensors have become an important tool to pre¬ 
cisely detect localized charges in nanoelectronic systems. 
A common technique to realize a charge detection scheme 
in two-dimensional electron gases (2DEGs) is placing 
a quantum point contact (QPC) close to the structure 
carrying charges. This technique has been successfully 
used to realize e.g. coherent spin and charra manipula¬ 
tion [Il,i or time resolved charge detection^, 1^. More¬ 
over, QPC based charge detection is a well-established 
technique to read out spin qubits in III/V heterostruc¬ 
ture 2DEGs il. It has further been shown that the 
detection fidelity can be enhanced by substituting the 
QPG for a single electron transistor (SET) which fea¬ 
tures a higher charge to no-charge signal ratio @ . 

The described detection schemes are even applicable to 
the relatively novel field of carbon based nanostructures, 
i.e. carbon nanotube ©[HI and graphene quantum dots 
(QDs) [Hill. Due to their weak spin-orbit [l3 . fl^ 
and hyperfine interaction [T^ . these materials promise 
long-living spin states making them potentially interest¬ 
ing for quantum computational applications fl^ . To 
detect the charge state of carbon nanostructures, nar¬ 
row graphene ribbons have been successfully integrated 
in carbon nanostructures and its functionality has been 
proven by a number of experiments [l^, [ 13 ) 113 ] • 

The electronic transport through graphene nanoribbons 
can be described by a stochastic Coulomb blockade, 
i.e. a series of QDs induced by edge and surface disor¬ 
der [ll, [l^ . Here, the slopes of the Coulomb resonances 
can be used to detect charging events in nearby quan¬ 
tum dots. This approach has been employed to per form 
charge sensing on individual graphene QDs [IlMliS- 
[^ . including time resolved detection [^ . 

While graphene QDs together with integrated charge de¬ 
tectors can be fabricated out of a single graphene sheet 


in a single lithography step, charge sensing in carbon 
nanotubes is more difficult to realize. In the past, the 
most common approach was to fabricate close-by metal¬ 
lic SETs [IjJSlor to couple a CNT based SET via a 
metal gate [^, . More recently, an all-carbon hybrid 

device has been demonstrated where charge states in a 
CNT have been detected by the current through a nearby 
graphene nanoribbon-based charge sensor 
Beside its ability to detect charge states, an equally im¬ 
portant property of a charge detector is its influence on 
the quantum dot itself. This so-called back action is cru¬ 
cial and needs to be understood if the detector is used 
in a system where a quantum dot has to be precisly con¬ 
trolled (e.g. in quantum computational schemes). 

Here, we report on the fabrication and characterization 
of graphene nanoribbon-based charge sensors integrated 
in graphene and carbon nanotube quantum dot devices. 
Effects of back action and counter flow are studied which 
give a clear indication of an energy transfer from the 
charge detectors to the quantum dots. 


II. DEVICE FABRICATION 

The fabrication of the device consisting of a graphene 
quantum dot with a nearby graphene charge detector is 
based on mechanical exfoliation of natural graphite. The 
flakes are deposited on highly doped Si’*''*' wafers with a 
290 nm thick Si02 layer. Single-layer graphene flakes are 
identified by Raman spectroscopy. Individual graphene 
flakes are subsequently nanostructured by electron beam 
lithography (EBL) and reactive ion etching (RIE) using 
an Ar/02 plasma. The devices are contacted by an addi¬ 
tional EBL step followed by metal evaporation. Fig.[T](a) 
shows a scanning force microscope (SFM) image of the 
investigated device. The graphene QD is located in the 
center of the image and has a diameter of around 130 nm. 
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FIG. 1. (color online) (a) Scanning force microscopy (SFM) 
image of a device consisting of a graphene quantum dot (QD) 
connected by graphene source (S) and drain (d) leads. Two 
charge detectors (CD) are realized by graphene nanoconstric¬ 
tions and positioned in the vicinity to the QD. (b) SFM im¬ 
age of a carbon nanotube (CNT) with a nearby graphene 
nanoconstriction acting as a CD. The bine areas indicate the 
position of the contacts to the CNT. The measurement cir¬ 
cuits are indicated by the schematic drawings in (a) and (b). 
(c) Condnctance of the graphene QD as a function of the back 
gate voltage Vbg in the Coulomb blockade regime showing 
distinct Coulomb resonances, (d) Measurement of Coulomb 
resonances of the CNT QD. 


The QD is connected to source (S) and drain (D) leads 
via narrow graphene constrictions acting as tunable tun¬ 
neling barriers. The charge detectors (CDs) are realized 
by two nanoribbons measuring a width of about 70 nm 
and positioned 50 nm from the QD. Apart from detect¬ 
ing charges the CDs can also be used as lateral gates if 
a voltage is applied with respect to the QD. Addition¬ 
ally, the underlying Si"'"’'" substrate can be operated as a 
global back gate tuning the overall Fermi energy of the 
device. 

Fig- mb) shows a SFM image of our second device scheme 
which consists of a CNT QD with a nearby graphene 
CD. To fabricate these CNT/graphene hybrid devices the 
CNTs are grown on Si+“'' /Si02 substrates by chemical 
vapor deposition (CVD) using a ferritin-based iron cat¬ 
alyst method [l^. This technique allows controlling the 
density of CNTs to be approx. 1 to 2 CNTs per The 
single-walled CNTs have a typical diameter of around 1.5 
to 2 nm and measure up to several micrometers in length. 
After the CNT growth, graphene is deposited on the same 
substrates by mechanical exfoliation. Graphene flakes in 
close vicinity to a CNT are selected and patterned into 
nanoribbons by EBL and RIE. Subsequently, they are 
contacted as described above. In the particular device 
investigated in this study (see Fig.[TKb)) the contacts are 
350 nm apart. Furthermore, the graphene constriction 
acting as the CD is located at a distance of 150 nm from 
the CNT and is 100 nm wide. 


III. MEASUREMENTS AND DISCUSSION 


Low temperature transport measurements are per¬ 
formed in a dilution refrigerator (at a temperature 
T < 100 mK) and a pumped '^He cryostat (T ss 1.5 K). 
Home-built low-noise DC amplifiers (amplification factor 
of 10®, bandwidth < 1 kHz) and low-frequency lock-in 
techniques are used to measure the current through 
the devices. Prior to the study of back action effects, 
suitable voltages are applied to the global back-gate and 
to the local lateral gates to tune the device into the 
Coulomb blockade regime. Please see Ref. [H, [13 for 
further details. Figs. [He) and 1(d) show the conductance 
of the graphene (c) and CNT (d) QD in dependence of 
the voltage applied to the back gate (Vbg). Distinct 
Coulomb peaks are visible which occur at every charging 
event of the corresponding QD. Fig. UKa) shows the 
current through the graphene QD as a function of the 
applied bias voltage {Vqd)- The observed diamond 
shaped areas of suppressed current constitute another 
distinct feature of the Coulomb blockade. Please note 
that the functionality of the nearby charge detectors of 
the particular devices shown in this study have been 
investigated in great detail and the results have been 
published in Refs. [H, 113) 113 • 

To study back action effects, we repeat the measure¬ 
ments of Coulomb peaks (cf. Figs. [TJe) and 1(d)) and 
bias spectroscopy (cf. Fig. [2Ka)) for a finite bias voltage 
applied to the CDs. Fig. [2{b) shows a bias spectroscopy 
measurement at Vcd = 20 mV in the same range as 
shown in Fig. [21 a). We observe an overall increase of 
the transmission through the QD for Vcd = 20 mV 
which is in agreement with earlier observations of back 
action effects in comparable devices [13 . Moreover, 
for some peaks the transition between positive and 
negative current moves away from the the line of 
Vqd = 0 mV (see dashed lines in Fig. |2{b)). For a 
detailed analysis Fig. He) shows corresponding line cuts 
at Vqd = 0, ±0.1, ±0.2, ±0.3, ±0.4 and ±0.5 mV. Here, 
the red trace illustrates the data set for Vqd = 0 mV. 
Evidently, the interaction with the CD (or back action 
of the CD) can result in current Iqd flowing against 
the biased direction which is most prominent for the left 
most peak. This effect is further enhanced at higher 
Vcd- a similar behavior can be observed for the CNT 
based device as shown in Eig. HKd). Here, the current 
through the CNT QD Iqd is shown for Vcd = 100 mV. 
Eor Vqd = 0 mV Iqd is negative and continuously 
develops into a positive peak as Vqd is increased to 1 mV. 


In order to discuss the observed behavior we employ a 
simple model of transport through quantum dots derived 
from rate equations considering only one ground state in 
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FIG. 2. (color online) (a) and (b) Current through the 
graphene QD Iqd as function of Vqd and Vbg mea¬ 
sured at VcD = 0 V and Vcd = 20 mV, respec¬ 

tively. (c) Line cuts through plot (b) at bias voltages of 
Vqd ~ 0, ±0.1, ±0.2, ±0.3, ±0.4 and ±0.5 mV. The zero at 
Vqd = 0 mV is highlighted in red. (d) Coulomb peak of 
the CNT QD in dependence of a voltage applied to the CD 
with respect to the QD Vq measured at Vbg = 46.2 V, 
Vcd = 100 mV and applied bias voltages Vqd = 0 mV 

to 1 mV as indicated in the panel. 


the QD: 

I = r[fs{E)-fD{E)] 

_ r r 

~ c‘(.Vo-Vg) + Vqd/2 c(Vo-Vg)-Vqd/2 ’ 

1 ± e 1 ± e '^BTD 

Here, F is the overall tunneling rate of the tunnel 
barriers and fs{E), fniE) are the Fermi distributions 
of the source and drain lead, respectively. Moreover, Vq 
is the center of the Coulomb peak, a the gate lever arm, 
Va the applied gate voltage and Vqd the bias voltage. 
The parameters Ts and Td describe the broadening of 
the Fermi functions fs{E), foiE) related to the two 
leads. A corresponding schematic of the model is shown 
in Fig. a). Here, the Fermi distributions fs{E) and 
foiE) are assumed to exhibit a different amount of 
broadening with Ts < Td and are offset by a positive 
bias Vqd- This directly results in two scenarios where 
the gate voltage can either tune the QD state to a 
chemical potential E where fs{E)-fD{E) > 0 (scenario 
1) or fs{E)-fD{E) < 0 (scenario 2). Scenario 1 thus 
results in an enhanced flow of electrons in bias direction, 
while in scenario 2 the electrons flow against the bias 
direction. 

Figs.EKb) and 3(c) show fits (red traces) of the described 
model to the experimental data (blue data points) at 
Vcd = 10 mV (b) and Vcd = 30 mV (c). Evidently, 
the model successfully reproduces the entire shape of 
the measured Coulomb resonances for different QD 
biases. In particular, the two scenarios 1 and 2, where 
the current is either enhanced (1) or flows against the 




FIG. 3. (color online) (a) Illustration of the model assuming 
an asymmetric broadening of the Fermi distributions fs{E) 
and fD{E) in the leads and only considering a single ground 
state contributing to electronic transport. The two red lines 
indicate possible positions of QD levels and the corresponding 
direction of net current is indicated by the green arrows, (b) 
and (c) Comparison of experimental data (blue) and fits (red) 
according to the model assuming single level transport. Data 
for bias voltages Vqd ranging from -0.4 mV to ±0.4 mV are 
shown in each panel. The bias applied to the charge detector 
Vcd measures 10 and 30 mV, respectively, (d) Extracted 
broadening parameters Ts (blue) and Td (red) as function 
of Vcd- Multiple data points per Vcd value correspond to 
different measurements in a bias regime IVqdI < 0.4 mV. 
The solid lines represent linear fits to the data, (e) Absolute 
peak current Ipeak extracted from the negative shoulder of 
each resonance (see arrow labeled by (1) in (b). The solid 
line represents a quadratic fit. 


direction of the applied bias (2) can be nicely explained 
(see e.g. arrows and labels in Fig. E^b)). To obtain 
more quantitative results, we fit the model to data sets 
measured for different Vcd and extract the parameters 
Ts and Td- Figure E^d) summarizes the results as a 
function of Vcd- We observe roughly a linear depen¬ 
dency of both effective temperatures Tg and Td on 
Vcd- However, while Td only increases from « 1.0 K at 
Vcd = 0 mV to « 2.5 K at Vcd = 50 mV, Tg shows 
a strong increase from Ril.OKtossSK within the 
same range. This finding clearly shows the increasing 
asymmetry of the Fermi distribution broadening of both 
leads. 

In addition, we analyzed the maximum current flowing 
through the QD. Fig. He) shows the absolute peak 
current extracted from the negative shoulder of each 
Coulomb resonance centered around Vc ~ 24.11 V as 
function of Vcd- A quadratic dependency of the peak 
current Ipeak on the detector bias is observed. This 
finding suggests that the current induced in the QD is 
approx, proportional to the electrical power dissipated 
in the CD. 

However, the model described above has a number of 
shortcomings. In particular it does not provide the ex- 
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pected electron temperature for zero CD bias, nor does 
it adequatly describe all observed transport phenomena. 
For example, Coulomb peaks exhibiting a total inver¬ 
sion of the current such as the one shown in Fig. (Sic) at 
Vbg ~ 24.08 V or Fig. [2][d) cannot be explained within 
the limits of the single level transport model. Thus, all 
these make it difficult to rely on and fully justify the use 
of this oversimplified model. 

Additional data, which cannot be described by our 
model is illustrated in Fig. lU^a). Here, we show a 
sequence of Coulomb peaks of the CNT QD recorded in 
a gate voltage {Vq) range of 500 mV for Vcd = 0 mV 
(black trace) and Vcd = 100 mV (red trace). Biasing 
the charge detector with Vcd = 100 mV leads to an 
enhancement of the current in regime A {Vc < 220 mV) 
by approximately one order of magnitude. In regime 
B [Vq > 220 mV) a reversal of the current occurs. 
Additionally, the Coulomb peak heights in both regimes 
are modulated by an envelope function with a period of 
roughly 300 mV. 

For the discussion of the observed behavior we employ 
a slightly more advanced model which considers three 
levels [i.e. the ground state (GS) and two excited states 
(ES)] contributing to transport through the QD. Similar 
to the single-level model above we further assume an 
asymmetric broadening of the Fermi distribution in the 
leads with effective temperatures Tg < Td and a finite 
applied bias voltage Vqd > 0- A schematic drawing 
of this model is shown in Fig. HKb). To describe a 
Coulomb resonance with completely inverted current 
flow, we furthermore assign different tunneling rates 
to the three different QD states. Considering higher 
rates for the two ES than for the GS (indicated by the 
widths of the arrows in Fig. 0Kb)) the net current flowing 
through the QD can be entirely negative as illustrated 
in Fig. UKc). Indeed, a comparison of the experimental 
data in Fig. [2Kd) to the model in Fig. 0Kc) shows an 
apparently qualitative agreement. 

Similar as for the single level model, values of Ts and 
Td for the leads of the CNT QD can be determined. 
Fig. HKd) shows the extracted parameters as function of 
Vcd- Ts and Td both increase approx, linearly with 
Vcd and show a qualitatively similar behavior compared 
to the values extracted for the graphene QD analyzed 
by the single-level model (cf. Fig.[3Kd)). The increase of 
Ts and Td with Vcd are measured to be ms « 49 K/V 
and Tod ~ 30 K/V, respectively. Thus they differ only 
by a factor of 1.6 in contrast to 6.1 for the graphene QD 
(cf. Fig.EKd)). 

If the tunnel rates of the three levels are assumed to 
be dependent on energy E, the three level model is 
furthermore capable to describe the sequence of Coulomb 
peaks shown in Fig. SKa). A corresponding sequence of 
calculated peaks for different broadening asymmetries 
{Ts = Td (black), Tg = 2 Td (blue) and Ts = 4Td) 
are shown in Fig. HKe). To obtain the illustrated traces, 
energy dependent tunneling rates fluctuating by over 
four orders of magnitude had to be assumed. The 
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FIG. 4. (color online) (a) Current through the CNT Iqd 
in dependence of Vg at Vqd = 0.3 mV, Vbg = 46.2 V, 
Vcd = 0 mV (black curve) and 100 mV (red curve). The 
behavior of the Coulomb peaks ranges from strong enhance¬ 
ment (regime A) to strong reversal (regime B) of the current, 
(b) Illustration of a model taking into account three states 
with different tunnel coupling (see different width of green 
and blue arrows) contributing to transport and an asymmet¬ 
ric broadening with Ts < Td- (c) Calculations according to 
the model illustrated in (b) for bias voltages ranging from 
0 to 1 mV. (d) Extracted broadening parameters Ts (blue) 
and Td (red) of the CNT QD as function of Vgd- The solid 
lines represent linear fits to the data, (e) Calculated current 
through a QD system with broadening parameters Ts and Td 
considering energy dependent tunneling rates which fluctuate 
by over 4 orders of magnitudes. 


calculated traces are in qualitative agreement with the 
experimental data [(compare Figs. 4(a) and 4(e)]. 

The single-level and three-level transport model dis¬ 
cussed above both rely on an asymmetric broadening of 
the source and drain leads as their most important as¬ 
sumption. The experimental evidence and the compari¬ 
son to the model presented in Figs.[3]and|3]show that this 
broadening has to be induced by the interaction of the 
GNR based charge sensors with the QDs. More insights 
can be gained from the quadratic dependency of the 
Coulomb peak current on Vcd (see Fig.[3Ke)) which sug¬ 
gests a proportionality to the power dissipated in the CD. 
For typical Vcd used in our experiments the dissipated 
power (which is most likely predominantly transferred to 
the CDs phonon bath) is on the order of hundreds of 
pW to tens of nW. However, the transfer of the thermal 
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energy from the graphene CDs to the QDs is hard to ex¬ 
plain. While detector generated phonons are regarded as 
a relevant origiir for back action in GaAs quantum dot 
devices [2^, phonon mediated energy transfer via the 
substrate is unlikely in carbon based devices due to the 
mismatch of the phonon spectra of graphene, CNTs and 
Si02. It is more likely that energy is transferred by pho¬ 
tons emitted by the current passing through the etched 
graphene nanoribbon (GNR) based detectors [1^. Pho¬ 
ton absorption in the leads of the QD devices could in 
turn lead to an effective broadening of its Fermi distri¬ 
bution. Moreover, it is possible that the occupation of 
the stochastic QDs formed in the GNR (sol - ls^ fluctu¬ 
ates. This would alter the potential landscape of the ca- 
pacitively coupled QD and allow charge pumping effects. 
Averaged over time this would have the same effect as 
an effective increase in temperature. Finally, the source- 
drain broadening asymmetry assumed in the model is 
likely to be a direct consequence of the asymmetry in the 
device geometry (see Figs. [Ha,b)) where one of the QD 
leads is closer to the GD than the other. 


IV. CONCLUSION 

Back action effects of graphene nanoribbon based 
charge detectors (CDs) on graphene and carbon nan¬ 
otube quantum dots (QDs) have been studied. We show 


that the current through the QDs is strongly influenced 
by the voltage applied to the CD. In both types of devices 
this voltage can lead to an enhanced transport through 
the QD but also to a reversal of the current. We employ a 
simple model considering transport only involving a sin¬ 
gle ground state which can explain the bias dependent 
behavior of isolated Coulomb peaks. A more advanced 
model taking into account excited state transport and 
state dependent tunneling rates proves suitable to fur¬ 
thermore describe longer sequences of Coulomb peaks. 
Importantly, both models assume an asymmetric broad¬ 
ening of the Fermi distribution in the leads of the QDs 
induced by the CD. Further experiments have to be con¬ 
ducted to obtain a better understanding of the mecha¬ 
nism of energy transfer between the charge detector and 
the QD. The reported study gives important information 
on the influence of graphene charge detectors on carbon 
based quantum dots which have to be considered when 
designing quantum dot circuits. 
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